Length dependent dynamics of microtubules 
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Certain regulatory proteins influence the polymerization dynamics of microtubules by inducing 
catastrophe with a rate that depends on the microtubule length. Using a discrete formulation, 
here we show that, for catastrophe rate proportional to the microtubule length, the steady-state 
probability distributions of length decay much faster with length than an exponential decay as seen 
in the absence of these proteins. 



INTRODUCTION 



Microtubules are important components of the cy- 
toskeleton of the cell These are hollow micro- 

scopic tubes formed by a, /3-tubulin hcterodimers. These 
dimers are joined end-to-end to form protofilaments with 
alternating a, /3 subunits. The wall of the microtubule 
has a staggered arrangement of typically 13 such protofil- 
aments. Experimental observations of the polymeriza- 
tion dynamics of microtubules reveal that microtubules 
alternate between persistent polymerizing and depoly- 
merizing phases. This unusual feature is known as the 
"dynamic instability" 0- The transition from a poly- 
merizing to a depolymerizing phase is known as "catas- 
trophe" and the same from a depolymerizing phase to a 
polymerizing phase is known as " rescue" . A popular idea 
regarding the cause of the dynamic instability is that it 
occurs as a result of two competing processes, the addi- 
tion of GTP bound tubulins at the growing tip and the 
hydrolysis of the GTP units of the tubulins present on 
the microtubule . The GTP bound tubulins provide 
a stabilizing cap to the growing end and if this cap is 
lost due to rapid hydrolysis of the GTP units, the micro- 
tubule enters into a depolymerizing phase. 

Microtubules are responsible for various intracellular 
organisation, exerting pushing or pulling forces etc. It is, 
therefore, important that the polymerization and depoly- 
merization of microtubules are appropriately regulated so 
that the length distribution of microtubules satisfies the 
requirement of the cell. It is found that the microtubule 
dynamics is linked to a large number of regulatory mech- 
anisms of microtubule associated proteins [l|. In vivo 
experiments show that a variety of proteins known as 
polymerases or depolymerases stabilize or destabilize the 
microtubules by reducing or enhancing the catastrophe 
rate @ respectively. Although the mechanism as exactly 
how the microtubule associated proteins influence the dy- 
namics is not known, there are experimental evidences 
of length dependent regulation of microtubule dynamics 
under different cellular contexts [^-Q. While in certain 
cases, there are speculations about the microscopic origin 
of such regulations Q , recent in vivo studies on two dif- 
ferent systems find depolymerases belonging to kinesin-5 
Q and kinesin-8 family of proteins responsible for 
such length-dependent regulation. In particular, live-cell 
imaging of interphase fission yeast cells @ suggests that 



kinesin-8 regulation leads to an increase of the catastro- 
phe frequency with the microtubule length. Since, in this 
case, the depolymerization is length dependent, cells may 
selectively control the (dc)polymerization dynamics of 
different populations of microtubules. As a consequence, 
the activity of such depolymerases is important for for- 
mation of ordered microtubule-based structures during 
cell division. Recent numerical simulations suggest that 
such length dependent regulation of polymerization dy- 
namics is crucial also for microtubule related processes 
during cell division [lo| . 

A possible mechanism through which kinesin-8 might 
introduce a length dependence in the catastrophe fre- 
quency has recently been proposed in ll|. The authors 
of [ll| simulate the regulatory mechanism based on the 
following postulates. The processive depolymerases, af- 
ter binding, walk towards the plus end of the microtubule 
and accumulate near its tip. Since these molecules can 
bind anywhere along the length with equal probability, 
and unbind rarely, it is likely that longer microtubules 
have larger density of the proteins at the tip than the 
shorter ones. Once the motor reaches the tip, it encoun- 
ters a 'push' from the subsequent motors. A 'pushed' 
motor at the tip finally breaks off taking the terminal 
tubulin along with. This gives rise to catastrophe due to 
the shortening of the GTP cap. These studies indicate 
that the catastrophe frequency increases linearly with the 
length of the the microtubule, as observed experimentally 
from live-cell imaging of fission yeast cells [9| . 

In the context of such dynamic instability, the cen- 
tral properties of interest are the length distributions of 
a set of microtubules [H, [13 ■ Apart from being useful 
for obtaining various statistical averages as a function of 
the rate parameters, these distributions lead to predic- 
tions of various quantities such as microtubule's mean 
life time, elongation time [l^ |. mean chromosome search 
time [l^ etc. that can be tested experimentally. In this 
work, by extending the two state model proposed earlier 
by Hill [13 , we study how the length dependent catastro- 
phe rate affects the length distribution of microtubules. 
Motivated by their own observation, authors of Q have 
proposed a length-dependent two-state model in which 
the catastrophe rate depends linearly on the microtubule 
length [l6| . Following a continuum approach, they obtain 
microtubule length distributions and use these distribu- 
tions further to obtain positional information regarding 
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the preferential cite for the future growth of the cell. In 
view of the relevance of the distributions, we are mo- 
tivated to adopt a discrete formulation to obtain exact 
expressions of various length distribution functions. The 
length distribution obtained in [l6j has a Gaussian form 
and our results agree with those of [l6| under a special 
condition satisfied by the growth and catastrophe rates 
and the microtubule length. 

The discrete description, unlike the continuum ap- 
proach involving two partial differential equations, re- 
quires two infinite sets of coupled first order dynami- 
cal equations. Since further consideration of the internal 
protofilament structure of the microtubule would lead to 
more complexity, we have considered a microtubule to be 
a polymer of a single protofilament as a first approxima- 
tion. Given the fact that the longitudinal bonds within 
a protofilament are strong in comparison with the lateral 
interaction [l3 , it is possible to extend the model by tak- 
ing neighboring protofilaments into account. However, it 
may be worth noting that considering a "gating rescue" 
effect from the neighboring protofilament as done for the 
composite model of leads to a behavior qualitatively 
similar to the single protofilament case. Other models of 
polymerization dynamics with len gth dependent attach- 
ment and detachment rates exist [l^, but these models 
do not take into account the persistent nature of growth 
and shrinkage. In view of this, the present analysis, al- 
lows a direct comparison with the results of [1^ in the 
steady-state. 

The evolution equations are written for and p~ 
which arc probability densities of having microtubules of 
length n in a growing or a shrinking phase respectively. 
The evolution equations for these probability densities 
are 



dt 



raPn-l + r-+Pn " + fa)p^ , n > 2, 



(1) 

= rdPn+1 + r+-np+ - {rd + r_+)p,T, n > 1. (2) 

Here, r^, rd, r-| and r |_ denote the tubulin attachment 

rate, tubulin detachment rate, catastrophe rate and the 
rescue rate respectively. As seen in these equations, the 

catastrophe rate, n, for transition from p'^ to p~ is 

proportional to the microtubule length, n. 



II. LENGTH INDEPENDENT CATASTROPHE 
RATE: CORRESPONDENCE WITH PREVIOUS 
RESULTS 

Here we consider the case where the catastrophe rate 
is independent of the length of the microtubule and find 
correspondences with the results of reference [l^. The 
steady-state recurrence relations satisfied by the proba- 
bility densities arc 

r-ari.i +r_+p^ - (r+_ +ra)p+ = for n>2 (3) 



^dPn+i + 'r+-pt - [Td + r^+)p„ = for n > 1. (4) 

Initially, we assume that p^ remains constant at a given 
value. We use the generating function method to find 
the dependence of the probability densities on the micro- 
tubule length. The generating functions defined as 



ri+2 



(5) 



Tl = 71=0 

satisfy the following equations 

[raX - (r+_ + ra)]G{x) + r^+H{x) = -rapt (6) 
[rd - {rd + r^+)x]H{x) + r+^xG{x) = (r^ -t- r_+)pj" 
^r+^pt- (7) 
Writing these equations in the matrix form as 

AG = B, (8) 

where 

Gix) \ 
H{x) j ' 



A = 



TaX - (r+_ + Ta) 



r-i X 



rd - {rd + 7-_+)a; 



and 



B 



-raPi 
{rd + r^+)p^ -r+^p+ 



we find the solutions for G{x) and H{x) as 

G{x) = r^I^i^rf " ("^d. + r_+)x}rapt + r+_r_+p+ 



-r-+{rd + r-+)p^ ], 

H{x) = —-^[rj^^raPiX + [VaX - {ra + ?+_)} X 



(9) 



{{rd + r^+)p^ -r+„_p+}]. 



(10) 



where | .A [= (1 - x)[ra{rd + r_+)a; - r<j(ra -f r+_)] is 
the determinant of matrix A. Clearly, matrix A is sin- 
gular at a; = 1 and this poses a question on the basic 
assumption of the steady-state. However, this singular- 
ity is a removable one if one considers the time evolution 
of the probability density, p^, and the probability den- 
sity of nucleating sites pq. With the following evolution 
equations [l^ . 



dpi 

dt 
dpo 

dt 



raPo + r-+Pi - {ra + r+-)p^ (11) 
-raPo + rdPi, (12) 



one obtains the steady-state conditions p^ = ^Po and 
p^ = DoPo, with Dq = I'Xll+lt] ■ It is straightforward 
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to see that the singularities in ^ and (|T0)) disappear 
once the ratio pf /Pi = r^^+r'' substituted in ([8|). 
Subsequent algebra leads to the steady-state probability 
densities 



Pn = -DoPo, and p„ = — Vo 



(13) 



as given in [12|. It is important to note that the ratio 
Pn/Pn as obtained here docs not match with the results 
of Appendix A of [l^. (Also see the discussion below.) 
The steady-state probability density p'^ decays exponen- 
tially with n as 



Pn 



Po exp[— an], with a~^hiDo, 



(14) 



if a is positive. In general, with constant p^, the singular- 
ity at X = 1 is removable provided Pi evolves with time 
such that the above ratio of pf /Pi is satisfied. This ratio 
is approached as t — oo with a correction disappearing 
exponentially with time. 



III. LENGTH DEPENDENT CATASTROPHE 
RATE 

In order to analyze the evolution equations with length 
dependent catastrophe rate, we consider the generating 
functions as given in ([5|) with time dependent probability 
densities. The generating functions are, in this case, time 
dependent. Summing ([1]) and ^ over n, we obtain the 
following differential equations for G{x,t) and H{x,t). 



dp^{t) dH{x,t) 

dt dt 
2r+^xG{x, t) + r+^x^G'{x, t) + 
{r+-pt - (rd + r^+)p^), and 



[rd-ird + r-+)x]H{x,t) + 



(15) 



dG{x, t) 
di 



raPt + {-ra - 2r+_ + rax)G{x,t) 



xG' {x, t) + r ^H{x, t). 



(16) 



Here, prime denotes derivative with respect to x. One 
could also directly consider the steady-state equations 
and obtain the probability densities by a transfer ma- 
trix approach. However, since its use is difficult in the 
length-dependent case, we follow the generating function 
method here. Using Eq. ^ for n = 1, all the x indepen- 
dent terms of equation (|15p can be together expressed in 
terms oi p^it). Taking Laplace transform of Eqs. (fTS)) 
and (froi). we have 



H{x, s)[sx - rd{l - hix)] = -rdP2{s) + 
2r+^xG{x, s) + r+^x^G'{x, s), and (17) 
G{x, s)[s + ra{l -x) + 2r+_] + r+^xG'{x, s) ~ 
r-+H{x,s) =rapt{s), (18) 



where hi = and /(x, s) = dt e-'* f{x,t). 

While obtaining these equations, we have assumed 
G{x, t = 0) = H{x, t ~ Q) — i.e. initially there are 
no growing or shrinking polymers of length n > 2. These 
two equations can be combined to obtain a single differ- 
ential equation for G{x, s) as 

G"(x, s)[(r+_s + r<jr+_)a;^ - r^r+^x] + rdr^+p^is) + 
G{x, s)[A + Bx + cx^] — ra[sx - rd{l ~ hix)]pi{s), 

(19) 

where 

A = -srd - raTd - 2r+_rrf, (20) 
B = + srd + sr-+ + ra-s + raVd + raTdhi + 
2r+_s + 2r+_rd, and (21) 
C = -ras - raTdhi. (22) 

In order to obtain the steady state behavior, we use the fi- 
nal value theorem which is based on the assumption that 
the steady-state generating function G{x) = G{x,t — > 
co) exists and is given by G{x, t —5- cxd) = lim,,^o sG{x, s). 
The same is assumed to be valid for other functions Pi{s) 
andp^(s). As s — >■ 0, Eq. pj]) becomes singular at a; = 1. 
However, as in the length independent case, this singular- 
ity is a removable one when the appropriate steady-state 
expressions for p^ and p^ are used. After these substi- 
tutions, the steady-state generating function satisfies the 
following equation 

dG 

{r+-rd)x— h (di - d2x)G(x) = da, where (23) 

ax 

di = Tdira + 2r+_), d2 = ra{rd + r_+) and 
d3^ptra{rd + r_+). (24) 

This inhomogeneous equation has a solution of the form 
G{x) = Gc{x) + Gp(x) having a complementary part 



Gc{x) = Cox '■+-''<i exp[- 



r+-rd 



(25) 



with Co as a constant and a particular integral for which 
a scries of the form in Eq. ([5]) can be considered. Eq. ([M)) 
has a singularity at x = and this results in a singularity 
of the complementary solution at x = 0. Since a catas- 
trophe rate proportional to the microtubule length is ex- 
pected to decrease the probability densities than those 
of the length independent case and a diverging solution 
for the generating function is unlikely, we choose Co = 0. 
The particular solution gives the probability densities as 



p+ = ds/di, 



^1 nm=i(di +mr+_rd) 
Pn+2 = :rPn+i for n > 0. 



(26) 

for n>2, (27) 
(28) 
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Eq. (|27|) . can be rewritten in the form 



+ ^3 



r(di/(r+_r,) + 1) 



di \{r+-rd)J r((n-l) + di/(r+_rd))'^^^^ 

where r(a;) is the usual T function of argument x. Hence, 
for large n, the exponential decay of ([T4| is modified as 

~ exp[— a'ji — nluTi], (30) 

with a' = In ^^-'^^ - 1. As Figs. [1] and [2] show, the de- 
cay of p+ of p~, dominated by an exp[— nlnn] type term, 
is much faster than that of the length independent dy- 
namics [l^. The ratio VnlVn-i^ however, remains same 
as that of the length independent case. In Figs. [3] and 
131 we have plotted the probability densities for different 
attachment rates. Fig. 0] shows an overall increase of 
p~ with the attachment rate. This seems to be a conse- 
quence of rapid polymerization and hence an enhanced 
transition of microtubules from growing to the shrinking 
phase. With the present form of distributions of growing 
and shrinking microtubules, the average length of micro- 
tubules is given hy n k, 26.4 for ra = 14.37 sec~^ and 
other parameter values same as provided in Fig. |31 

In a continuum description of the microtubule length, 
one may have the following steady-state differential equa- 
tions [iq 

-raa%^ + r_+p- (L) - r+_ -p+ (L) = 0(31) 
oL a 

rcia%^+r+--p+{L)-r^+p-{L)=Q, (32) 
oL a 

where a is the length of the tubulin dimcr. In case of 
bounded solutions for p'^{L) and p^{L), it is found that 

^-Tfj^ — [3 ■ J^l- (|28p and also the corresponding cq. 
in the length independent case ( Eq. (fT3|) ) arc different 
from this. Although ([25)1 and its continuum analogue 
mentioned above are simple balance equations in terms of 
the rate parameters, the subtle difference between them 
becomes important here unlike the length independent 
case. In order to find the reason for such difference, wc 
do a Taylor expansion of Eq. (j28p in the dimcr length, a. 
Retaining terms up to first order in a, the ratio appears 



Replacing the first 



p+(L) - rA-^ "p+(i) dL 

order term using Eq. (I3ip. it can be seen that the ratio 
becomes same as that of [la| provided r+_/ra<<a/L = 
1/n. Rather than referring this as a condition for validity 
of the continuum approximation, it is more appropriate 
to mention that this condition expresses a competition 
between the catastrophe and the growth rate which, as we 
shall see below, distinguish two regimes with two different 
kinds of probability distributions. 

That the distribution in ^7} finally becomes a Gaus- 
sian distribution for /r^ << a/L, can be seen by 

approximating the logarithm of the product in the de- 
nominator of (|27p in the following way. 

n-2 

In ]^ {di + mr+^rd) = 




100 120 140 160 



FIG. 1: Steady-state probability densities for growing micro- 
tubules are plotted with length. Solid and dashed lines cor- 
respond to length independent and length dependent catas- 
trophe rates respectively. Parameter values used for the 
plot are — 10 s^^, rd = 40 s"^, r_+ = 0.02 s^^ and 
r+_ = 0.19 s-^. 
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FIG. 2: Steady-state probability densities for shrinking mi- 
crotubules are plotted with length. Solid and dashed lines 
have the same meaning as in Fig [1] Parameter values are 
same as those of Fig. [T] 



n-2 



(71 - 2)ln(ii + ^ ln(l 



771 — 1 
L 



(n - 2)ln(ii + / dx ln(l 

"'0 



{n - 2)lnrfi + 



{ra + 2r+_)a Jo 



I— « 

di 

a di 

L 

dx X. (33) 



Since /ra « a/L, we have used ln(l + x) sa x to 

approximate the integrand. It is straightforward now 
to see that ([55)1 finally leads to a Gaussian distribu- 
tion for the length. exp[— nlnn] form of Eq. (j30|) . on 
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the other hand, follows in the opposite situation where 

\/n = a/L « r_| /r^. This analysis shows that the 

higher is the growth rate the quicker is the fall of the 
probability distribution of growing polymers. This seems 
to be a consequence of the fact that a rapid growth in- 
duces more frequent catastrophe. 



FIG. 3: Steady-state probability densities for growing mi- 
crotubules are plotted as a function of length for different 
attachment rates. Values of parameters other than ra are 
rd = 37.5 r+- = 0.014 s"^ and r-+ = 0.044 
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IV. SUMMARY 



It has been found that the catastrophe frequency 
becomes dependent on the microtubule length due to 
the activity of depolymerases belonging kinesin-8 fam- 
ily. Here, we consider a two-state model of growing and 
shrinking microtubules having a catastrophe rate that in- 
creases linearly with the length of the microtubule. Fol- 
lowing a discrete approach, we study how the steady state 
probability densities and p~ of growing and shrink- 
ing polymers, respectively, of length n, vary with length. 
Our results show that, for large n, the decay of the prob- 
ability densities with n is dominated by an exp[— nlnn] 
type of term as compared to an exp[— an] kind of decay 
in the length independent catastrophe case. 



FIG. 4: Steady-state probability densities for shrinking mi- 
crotubules are plotted with length for different attachment 
rates. Parameter values are same as those of Fig. [S] 
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